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ABSTRACT 

In recent years it has been shown that the tidal coupling between extrasolar planets and their stars 
could be an important mechanism leading to orbital evolution. Both the tides the planet raises on the 
star and vice versa are important and dissipation efhciencies ranging over four orders of magnitude 
are being used. In addition, the discovery of extrasolar planets extremely close to their stars has 
made it clear that the estimates of the tidal quality factor, Q, of the stars based on Jupiter and 
its satellite system and on main sequence binary star observations are too low, resulting in lifetimes 
for the closest planets orders of magnitude smaller than their age. We argue that those estimates 
of the tidal dissipation efhciency are not applicable for stars with spin periods much longer than 
the extrasolar planets' orbital period. We address the problem by applying our own values for the 
dissipation efficiency of tides, based on our numerical simulations of externally perturbed volumes of 
stellar-like convection. The range of dissipation we find for main-sequence stars corresponds to stellar 
Q* of 10^ to 3x10®. The derived orbit lifetimes are comparable to, or much longer than the ages of 
the observed extrasolar planetary systems. The predicted orbital decay transit timing variations due 
to the tidal coupling are below the rate of ms/yr for currently known systems, but within reach of an 
extended Kepler mission provided such objects are found in its field. 

Subject headings: Convection; Turbulence; Methods: numerical; Planets and satellites: dynamical 
evolution and stability; Planet-star interactions; planetary systems 



1. INTRODUCTION 

With the discovery of several giant extrasolar planets 
that are extremely close to their stars, it has become clear 
that the dissipation of the tides that the planet raises on 
the star plays an im portant role in determinin g the or- 



bital evolution (c.f. Dobbs-Dixon et al . 
Burrows||2009l |Pont||2UU91 [Jackson eTa 



Matsumura et al. 2008 an iong others 



2004 



Ibgui fc 



already has a history (e.g . 



et al. 1998 Sasselov 2003 



Rasio et al. 1996 



The problem 



Ogilvie & Lin 



3961 iT erquem 
pM| , which 



led to the basic conclusion that a comprehensive theo- 
retical understanding of turbulent dissipation in stellar 
convection zones is lacking. The recent discovery of the 
the very close- in giantplanet WASP-19b made that very 
clear (Hebb et aV'20To\- 

In some cases the dissipation of the tides on the planet 
is important as well, but (at least for single planet sys- 
tems) this eventually leads to the planet being on a cir- 
cular orbit with its spin synchronized with the star, at 
which point the tides on the planet become independent 
of time and hence no dissipation occurs. The tides raised 
on the star will also have the tendency to synchronize the 
rotation of the star with the orbit, but in the case of plan- 
ets usually the angular momentum of the planetary orbit 
is not enough to achieve that. 

Commonly the rate at which energy is dissipated from 
the stellar tides is parametrized in terms of a quality fac- 
tor (Q*) defined as the fraction of the tidal energy dis- 
sipated in one tidal period. It is then assumed that this 
value depends only on the star, and not on the frequency 
of the tides and its value derived from observations of 



main sequ ence solar mass binary stars ( Meibom & Math- 
ieu||2005|) is used: 10^ < Q* < lO"^. This dissipation. 



however, is found to lead to unreasonably short lifetimes 
for at least two systems — Wa .sp 18 b ( jHellier et al.^2009| 
and WASP-19b ( |Hebb et alT]|2010[ ). i-br these systems 
the time it would take tor the planets to plunge into the 
star, if the above values are assumed, is about three 
orders of magnitude smaller than the estimated ages of 
the systems and four orders of magnitude smaller than 
the total lifetime of the star, making it extremely un- 
likely that such a sys tem should ever have been observed 
(Hellier et al. 2009), unless one assumes a mechanism 
that continuously resupplies these extremely close to the 
star regions with fresh planets. 

There is however a fundamental difference between the 
tides in a binary star and a star-planet system, namely 
that the mass of the secondary in the case of stars is large 
enough to spin the primary to synchronous rotation with 
the orbit, and in the case of planets it is not. And in- 
deed for all transiting planetary systems for which the 
spin period of the star (p^) is known it is found to be 
much longer than the orbital period (P). This distinc- 
tion is important, because when |P| > |p*/2| the time 
variable tidal perturbation can resonantly excite inertial 
waves in the star, thus resulting in several orders of mag- 
nitude larger shear, and hence dissipatio n, compared to 



the static tid e (Savonijc & Papaloizou 1997 Dintrans fc: 
|Ou yed|200 H fOjilvie fc' Lin.2004,.Wu 2005aIbl|Papaloizou 



& lyanov 2005; Ogil vie fc: Lin|2007 Ivanov fc Papaloizou 



pO?) iOgilvie||2009| ) 

Fot stars with surface convective zones the turbulent 
convective motions are thought to be causing a cascade 
of energy to occur from the large scales of the tides to 
smaller and smaller scales until eventually the finite vis- 
cosity of the plasma becomes important and converts this 
energy into heat. The usual treatment of this compli- 



2 



cated process is to assume that it b ehaves like an effective 
turbulent viscosity coefficient (c.f. Zahn||1966 Goldreich 



&: Nicholson||1977[ ). However, the different prescriptions 
have been diificult to r econcile with each othe r and with 

different observations (Goodman & Oh f997). 

( ,2009a) ; Pcncv et al., (2UU8 , 2009b 



mass and radius (see Sasselov||2003 for more details) 



Penev et al. 



20071 



used numerical simulations of stellar-like convection to 
show that the dissipative properties of the turbulent con- 
vective flow are indeed well approximated by an effective 
viscosity coefficient. Therefore the simulations allow us 
to derive its theoretical value. In this paper we combine 
these new results to construct a complete prescription for 
the effective viscosity and show that the resulting dissi- 
pation, in the absence of resonantly excited tidal waves, 
produces much higher values than those found obser- 
vationally for main sequence stars, and that those values 
are consistent with even the strongest current constraints 
on the tidal dissipation efficiency. In addition we derive 
the magnitude of the expected transit timing variations 
(TTV) due to orbital decay and compare them to possi- 
ble observational contraints (e.g. by the Kepler mission) 
that could provide a direct test of the proposed viscosity 
prescription. 

The organization of the paper is as follows: in Section 
[2] we introduce what is usually understood by the 
parameter, in Section [3] we discuss how we arrive at the 
turbulent viscosity we use to calculate orbita l decay, in 
Section [4] we follow |Scharleman"ii| ( |1981[ |1982|) to convert 
our turbulent viscosity to a tidal torque, m Section [5] we 
present the stellar structure and evolution models used 
in the calculation of the torque, in Section |6] we calculate 
the effective Q* which corresponds to this torque, the 
TTVs that our estimate of the dissipation would result 
in and the future lifetimes of close in extrasolar planets, 
finally in Section [7] we summarize our results. 

2. CLASSICAL APPROACH TO Q 

The tide raised by a planet on its star is a quadrupole 
wave of amplitude h; the tidal motions are of order 
(uj — fl)h which set up shear inside the star. The or- 
bital angular frequency, uj, exceeds the stellar spin one, 
i7, in the case we consider here, and the tidal forcing 
due to u! occurs with timescales that are much shorter 
than the correlation time of the turbulence in the star's 
convection zone. 

Dissipation in the star causes the tidal bulge h to lag 
behind by an angle 6, which is determined by the amount 
of coupling between the tide and the source of the dissipa- 
tion, presumably the turbulent eddies in the convection 
zone. One could compare the response of the star to that 
of a forced harmonic oscillator and relate 6 to a specific 
dissipation funct ion Q = {oj — fl)Eo/E — 1/2S (Murray 



fc Dermott|l999 ). This is how the tidal dissipation qual- 
ity factor (J* is defined, with Eq being the energy stored 
in the tidal bulge, and E, the rate of viscous dissipation 
of energy. This is another way to determine the lag angle 
S. 

The turbulent viscosity is introduced in the sense of 
Rayleigh-Benard incompressible convection, i't = \vl ~ 
/^/r, for a low forcing frequency. Here, I is the convec- 
tive mixing length, v is the convective velocity, and r is 
the convective timescale (eddy turnover time) . Then ba- 
sically Qs, oc GM/r^^uivt, where M and r* are the stellar 



3. EFFECTIVE TURBULENT VISCOSITY 

We need to c ombine the perturb atively derived effec- 
tive viscosity of Penev et al.| ( |200 9b) base d on realistic 
low mass star convective models with the IPenev et al.l 
(|2009a) direct viscosity from simulations with external 
forcing in order to get a complete and reliable prescrip- 
tion for the full v iscosity tenso r. 



As discussed in Penev et al. ( 2009b I the viscosity ten- 
sor is specified completely from five independent quanti- 
ties: A'o, Kq', Kqq', Ki, K -^. Two of th e se (K ^ and K2) 
were calculate d directly in 



the other three Penev et al. 



Penev et al. (2009a), and for 



] 2009 b ) provme perturbative 



values. In terms of these quantities the time and volume 
averaged rate of energy dissipation d ue to the turbulent 
flow is given by ( Penev et al.||2009b equation 13): 



^uisc (^) 



dz [ KiAl 4- K2AI + KoAl + 
+Ko,Al, +2K(,o,AaAa'], (1) 



where the |A„ip quantities are root mean square shear 
components defined by: 
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(2) 
(3) 
(4) 
(5) 



with the angle brackets denoting a time average. 

Equation [T] is simply the volume and time average of 
the rate of work done by an anisotropic viscous force on a 
stratified fluid subject to some externally imposed shear. 
The particular form of Eq. [l] assumes that locally the 
viscosity tensor is invariant under rotations around the 
axis of gravity (z), which must be true if stellar rotation 
is ignored. 



The directly obtained effective viscosity of Penev et al. 
(|2009a| and the perturbative estimates based on realis- 



tic low mass star simulations (Penev ct al. 2009b| ) and on 
an idealized simulation (Penev^ 3t al. 2008) all produce a 
linear scaling of the effective viscosity with period, for 
the range of periods available to those simulations. As 
is discussed in these works the linear scaling is not ex- 
pected to hold for arbitrarily small periods, because at 
such timescales Kolmogorov cascade should be a good 
approximatio n to the flow and in that case IGoodman] 
& Oh (1997) show that the loss of efficiency should be 
quadr atic with period. I n addition, the direct calcula- 
tions ( Penev et al.||2009a ) show that the effective viscos- 
ity saturates at long periods. 
To accommodate these three scalings, we will assume 
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the following form for the viscosity coefBcients 

2 

s 



K„ 



: mm 



T I T 



1 



2\l/2 ; 



with: 



n 



(6) 



(7) 



where, I is the mixing length, and IVmin and Umax are the 
dimensionless periods between which the linear scaling 
applies, Sm are the set of slopes with period of each effec- 
tive viscosity component, and are the corresponding 
zero crossings. 

This expression has been chosen to be continuous with 
period (P), reproduce the quadratic scaling we expect 
to occur at small periods, the linear scaling at periods 
comparable to the local convective turnover time (r) with 
a possible offset (A^), and a constant effective viscosity 
for muc h longer perfods. 



From Penev et al. (2009a): 



si = 0.084, S2 = 0.055. 



(8) 



Assuming that the scaling betwee n the slopes of these 
components and the rest from the Penev et al. (2009b I 
perturbative calculation holds, we arrive at the toUowmg 
values for the remaining three s values: 



SO' = 0.23, So = 0.07, soo' = 0.02 



(9) 



The offsets, A™, are less well constrained from the 



Penev et al. (2009a) and Penev et al. (2009b) simula- 
tions, because on one hand their values from the direct 
calculation depend on the method used to extract the 
effective viscosity and on the other the perturbative cal- 
culation does not predict any offsets so it is not clear of 
how to find them for the remaining components. 

For the components not corresponding to radial shear 
we will assume no offset. The Ai component is expected 
to be between the values d erived by the two direct meth- 
ods of Penev et al. (2009a) for a forcing strength of zero. 
Ignoring the ditterence between the weak forcing and a 
zero forcing c ase, which according to figure 13 of [Penev] 



et al. (2009a I is likely to be small, we conclude that: 



0.023 < Al < 0.061. 



(10) 



For Ao' we will show results with two different assump- 
tions: 

Ag, =0, Aj, =2.6Ai, (11) 

The last case corresponds to assuming that the offset 
scales the same way as the slope between the m = 1 and 
m = 0' components. 
Finally we need to specify Umax and Tlmin- From 



Penev et al. (2009a) we see that Umax — 2.4 is reason- 
able. i'he value of ^\^min is not seen in the simulations, 
but clearly Timiri < 0.3, so we will consider two cases 

Hmm = 0.1 and limin = 0.01. 

4. TIDAL TORQUE 

An analytical expression for the tidal velocity and the 
associated tidal torque w as derived with the sm allest 
number of assumptions by Scharlemann ( 1981 1982^ for 



a circular orbit and close to synchronous internal rota- 
tion. Here we will include the relevant results and refer 
t o those works f or thei r derivation. 

Scharlemann (1982) equation (17) gives the tidal 



torque as: 



T 



967r 



(12) 



where, is the radius of the star, is the maximum de- 
viation of the equipotential surfaces from spheres, n is the 
orbital angular frequency, Ai can be expressed through 
integrals of the viscous force over the star, ao — aoc is 
the departure from synchronous rotation, and u)2* and 
ao describe the differential rotation assumed to have the 
form: 



r2(r, 6) = nuj2* ao + 



I 



-x^ sin^ ( 



(13) 



The second term in the above expression does not enter 
in Al and so will not play an y further role. 



The value of Ai is given by Scharlemann ( 1982 ) equa- 
tion (16): 



Al 



20 



'92 



25 
28 



^04 ^ 



14 



(14) 



where the indices, 2 and 4, refer to the P|(cos0) and 
P|(cos0) associated Legendre polynomial components 
respectively of the following integral s of the viscous force 
(P) acting on the tidal velocity field ( Scharlemann ( 1982 1 
equations (5 c-e)): 



-dr. 



le{0) -- 



9 

■tanflPfl 



dr, 



' sin 6iPrf, 
9 



-dr. 



(15) 
(16) 
(17) 



5. STELLAR MODELS 



The viscosity tensor enters in the calculation of Ai only 
through its 8* order moments: 



'Km{P,r)dr, 



(18) 



The depth dependence of Km arises from the fact that 
we expect Km to scale as plvjT) (p, / and v are the den- 
sity, the mixing length and the convective velocity re- 
spectively) , all of which are functions of depth and it also 
depends of the convective turnover time, approximated 
as //f, which also varies with depth. 

In order to evaluate the right hand side of equation [18] 
we need p, I and v in the interior of the star. We will take 
these quantities from two stellar models: a 1.4M0 and a 
O.8M0 main sequence solar met alicity models calcul ated 
using the code described in Cody & Sasselov ( 2002 ) . 

The reason for choosing these two models is that 
presently almost all the known transiting extrasolar plan- 
ets are around stars that fall within this mass range, so 
our results will span the range of currently known sys- 
tems. 



4 




P[ms/yr] 



Fig. 1. — The tidal quality factor corresponding to the combined 
viscosity from Pcnev et al. (2009a); Penov et al. (2009b, 2008) for 
the range of assumptions described in Section |3] i'he blue dashed 
curves correspond to the IAMq model, the remaining curves are 
for the O.SMq model: Tlmin = 0.1 — green, Umin = 0.01 — red 



6. ORBITAL EVOLUTION 

From the tidal torque we can calculate the rate at 
which the orbit of a circularized and synchronized planet 
orbiting a slowly rotating star shrinks. For comparison 
purposes we converted this rate to an effective value, 
shown in figure [l] 

As we can see, for both stellar models and all as- 
sumptions of Section [31 the dissipation efficiency is much 
smaller than the usually assumed values of 10^ < < 
10^, implied observationally for main sequence binary 
stars. As discussed in the introduction we expect that 
this is due to the fact that in the case of planets the tides 
have a frequency that is too high to resonantly excite in- 
ertial waves in the star. 

The reason for twice as many curves being visible for 
the O.8Af0 case than for the 1.4M0 is that the two dif- 
ferent assumptions for Tlmin are indistinguishable for the 
higher mass model due to the fact that the convection in 
the high mass star is much more vigorous and the orbital 
period does not get below 0.1 convective turnover times 
for most of the convective zone for the range of separa- 
tions considered. The dissipation being more efficient for 
the low mass star is mostly due to the fact that it is much 
less centrally concentrated and hence a lot more mass is 
subject to the tides than for the high mass case. 

In figure [2] we show the TTVs that will be produced by 
this turbulent viscosity. Each colored region corresponds 
to the range of viscosity assumptions listed in Section 
[3) blue for the 1.4M0 model, and red for the O.SMq 
model and the red to the most dissipative. For reference 
we have added the currently known transiting extraso- 
lar planets that fall within the plotted range, marking 
the ones most affected by tides separately: WASP-12b 
- square, WASP-18b - rhomb, OGLE-TR-56b - circle, 
WASP-19b - up triangle, WASP-33b - down triangle. 
As can be seen, for even the closest and most massive 
extrasolar planets the rate of period change is less than 
a ms/yr. 

The small TTVs of course are to be expected since 
planets which decrease their orbital period quickly will 
exist only for a short amount of time (relative to the 
main sequence life of their parent star) before falling in- 




Fig. 2. — The rate of change of the orbital period of the planet in 
ms/yr. The colored regions correspond to the range of assumptions 
discussed in Section[3[ blue for the 1.4AfQ model; red for the O.SMq 
model. The black circles and the green symbols correspond to 
the currently known transiting extrasolar planets that fall in the 
plotted range. For each stellar model, the values of the TTVs for 
the regions going from left to right are 10, 1, 0.1, 0.01 and 0.001 
respectively. 



side the Roche lobe and being tidally destroyed, making 
the probability of finding such planets small. They are an 
order of magnitude smaller than what was predicted for 
OGLE-T R-56b using th e older linear prescription for dis- 
sipation ( Sasselov||2003 1 . However, if such close-in giant 
planets are discovered by Kepler, an extended multi-year 
mission could achive sufficient precision to constrain the 
theoretical work used here. 

To address the question of tidal destruction of planets 
we used the expression for the tidal torque to calculate 
the coupled evolution of the planetary orbit and the stel- 
lar spin for systems that are assumed to start with a 
non-spinning star and a planet in a circular orbit. We 
assume that the planet is synchronized, so we do not 
consider the tides on the planet and evolved the system 
forward in time until the dissipation drives the planet 
close enough to its star to be tidally destroyed, thus get- 
ting an estimate of the lifetime of the planet during which 
it will be observable. 

These lifetimes for a range of planet masses and sep- 
arations are shown in figure [3] for the two stellar mod- 
els discussed above. Each pair of blue curves shows the 
most and least dissipative effective viscosity assumptions 
for the 1.4M0 model and the red curves show the range 
for the 1.4M0 model. The lines to which the curves 
for a given model asymptote in the upper right corner 
correspond to the critical planet-star separation beyond 
which the planetary orbit has sufficient angular momen- 
tum to eventually spin up the star to synchronous rota- 
tion thus halting tidal evolution. Again for reference we 
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Fig. 3. — The future lifetimes of single planet systems on a cir- 
cular orbit around an initially non-rotating star for the two stellar 
models considered in the text: O.SA/q — red, 1.4Mq — blue. 
Each colored region corresponds to the range of assumptions from 
Section |3] The black points and green symbols correspond to the 
currently known transiting solar planets that fall within the plot- 
ted range of masses and separations. The remaining lifetimes for 
each model corresponding to the shown regions are 0.01, 0.1, 1, 10 
and 100 Gyrs going from left to right. 



have shown the known transiting planets. 

As we can see most of the currently observed planets 
have tidal lifetimes of many billions of years, and (af- 
ter considering Wasp-12b, Wasp-18b, Ogle-TR-56b and 
WASP-19b individually) we see that no planet has a life- 
time of less than a few hundred million years, the short- 
est lived being Wasp-18b with a lifetime ranging between 
100 and 1 60 Myrs for a 1.25 Mf:^ stellar model (the mass 
quoted in [Hellier et aL]|2009[ ) . 



7. DISCUSSION 

In this paper we have presented a direct theoretical 
calculation of the dissipation efficiency of tides raised on 
stars by a short period giant planet. The estimates of 
this efficiency are based on direct simulations of the dis- 
sipation of externally driven perturbations in a small, but 
still significantly stratified, piece of a convective zone. 

While due to numerical limitations we are unable to 
explore the entire range of the ratio of tidal frequency 
to local convective turnover time that occurs in the sur- 
face convective zones of low mass stars, we show that 



the rate of dissipation of tidal energy averaged over the 
entire convective zone is constrained to within a factor 
of ten. This is a significant improvement over the cur- 
rent range of three to five orders magnitude that many 
authors wishing to calculate the tidal ev olution of ex- 
oplanet orb its are forced to consider (e.g. iHellier et al. 
2009| [Jackson et al. 2009; Barker & Ogilvic 2009T|Iev" 



rard et al.|[2009i ,Patzold et al.,,2004, .Adams fc Laughlin 



2006F 



Since the usual way to parametrize the dissipation of 
tides is using the stellar quality factor Q*, we calculate 
effective values for this parameter as a function of the 
orbital frequency. We find that our estimates lie towards 
the upper end (small dissipation) of the usually assumed 
range for this parameter -10^<Q*<10^~ f^nd are 
significantly larger than the value needed to explain the 
fact that the maximum period to which solar type binary 
stellar orbits are circularized increases noticeably over 
the main sequence lifetime of these systems. 

We argue, however, that this is not necessarily a con- 
tradiction, since in tha case of binary stars synchroniza- 
tion of the stellar rotation to the orbit happens much 
faster than the circularization of the orbit. This means 
that the tidal frequency is close to the rotation rate of 
the star, for most of the time circularization occurs. This 
opens the possibility that inertial waves are resonantly 
excited in the stars, which results in a much enhance dis- 



sipation ([Savoni je fc P apaloizou 1997; Dintrans fcOuyed 
200T1 [ Ogilvie'fc Lin '20^4; Wu 2005a b; Papaloi zou H 

Ivanov 



Ivanovi.2005 
20071 lOgilviel 



Papaloizou 



Ogilyie fc Lin„2007; , 

|2009p , over what is likely to occur in planet 
which 



star systems, which typically have the star rotating much 
slower than the tides raised by the planet. 

This idea seems to be confirmed by a number of tran- 
siting extrasolar planets which are found to lie very close 
to their host star. So close in fact that were the value of 
implied by main sequence binary stars applicable to 
star-planet systems, we would conclude that we caught 
several of those pla nets in the last less than .1% of their 
lifetime ( Hellier et al.|2009[[Hebb et al.|2010l e.g.), which 
seems rather unlikely given the current sample of about 
a hundred known transiting planets. 

We show that our much smaller dissipation results in 
future lifetimes for all the currently known transiting 
planets comparable to their age (Fig. [s]), and is hence 
consistent with those observations. 

We also calculate the inferred TTVs due to tidal decay 
of the orbit (Fig. [2]) and find values of < Ims/yr for all 
currently know transiting planets. Much smaller than we 
can hope to observe from the ground in the near future, 
but perhaps within reach of an extended Kepler mission if 
a very short period massive planet happens to be among 
its targets, leaving the possibility to measure directly. 
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